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Abstract 

By measuring Drell-Yan massive lepton-pair's angular distributions, we can 
identify the polarization of the virtual photon of invariant mass Q which 
decays immediately into the lepton-pair. In terms of a modified QCD fac- 
torization formula for Drell-Yan process, which is valid even if Qt S> Q, we 
calculate the massive lepton-pair's angular distributions at large Qt- We find 
that the virtual photons produced at high Qt are more likely to be trans- 
versely polarized. We discuss the implications of this finding to the J/V' 
mesons' polarization measured recently at Fermilab. 

PACS Numbers: 12.38.Bx, 12.38.cy, 13.85. Qk, 14.70.Bh 
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Recent data on J/V' polarization measured by CDF collaboration at Fermilab Tevatron 
seem to be inconsistent with the predictions from various models oi ^/%p production [Q. The 
CDF data shows that the J/V' mesons produced at = 1.8 TeV become more longitu- 
dinally polarized as the transverse momentum {Qt) increases [|l|. On the other hand, the 
various theoretical calculations predict the J/V' mesons to be more transversely polarized at 
large Qt §• 

In order to help solving this puzzle, it was proposed to measure the Drell-Yan massive 
lepton-pair's angular distributions at the kinematics similar to that of J/?/' production, and 
to extract the virtual photon's polarization from the lepton-pair's angular distributions [^Q. 
Production of the virtual photons in the Drell-Yan process at large Qt and small has a 
lot in common with the production oi ^/ip mesons at high Qt- They both have two large 
physical scales: Qt and Q^, which is equal to Mjy^ in the case of i /ip production; and Qt 

is much larger than Q^ . If the collision energy is large enough, and the logarithmic 
contributions from resumming all powers of \n"'{Q\/Q'^) dominate the production cross 
sections, the virtual photon and J/V' production will share the same partonic subprocesses, 
labeled by the R in Fig. |l|. Then, the only difference between the virtual photon and J/V' 
production at large and high Qt is the difference in the fragmentation functions. The 
virtual photon fragmentation functions are completely perturbative if Q^ 3> ^'qcy) [liil, 
while the parton to J/?/' fragmentation functions involve final-state nonperturbative soft 
interactions 0. The measurements of the virtual photon and J/?/' polarization at high Qt 
can help us to isolate the final-state effect in J/-?/' production, and to narrow the questions 
about ^/ip formation. 

In this letter, we calculate Drell-Yan massive lepton-pair's angular distributions in Quan- 
tum Chromodynamics (QCD). We present the polarization parameter ody as a function of 
Qti which can be in principle measured at Fermilab Tevatron. We find that the virtual 
photons at high Qt are likely to be transversely polarized, which is different from the data 
on the polarization parameter aj/^ measured at Fermilab If our predictions for the ody 
are confirmed by future experiments, we may conclude that the non-perturbative final-state 
interactions are very important for the formation oi ^/ip mesons. 

When Qt and Q^ are both large, cross sections for Drell-Yan massive lepton-pair pro- 
duction in hadronic colhsions, A{Pa) + B{Pb) 7*(^ ^KQ)) + ^5 can be systematically 
calculated according to QCD factorization theorem Since we are mainly interested in 

the cross sections at high Qt and low Q^, we neglect the Z channel contributions in the 
following discussions. 

After integrating over the azimuthal angle distribution, the Drell-Yan lepton-pair's an- 
gular distributions, as sketched in Fig. ||, can be expressed as @] 



(^'^AB^U{Q)+X _ AB^U{Q)+X 



dQ"^ (IQt dy d COS 6 dQ'^dQj.dy 2(3 + ody) 



1 + ody cos 9 (1) 



where the polarization parameter ody is given in terms of the cross sections for producing 
a polarized virtual photon, 

d(TAB-.'rl{Q)+X ^^ j I d(TAB^-,*^{Q)+X _^ daAB^'yl{Q)+x '\^ 




dQTdy dQTdy j j \ dQ^dy dQrdy j 
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The subscripts, T and L in Eq. (||) represent the transverse and longitudinal polarizations, 
respectively. The inclusive Drell-Yan cross section in Eq. (P can also be expressed in terms 
of the cross sections for producing a polarized virtual photon. 



da 



AB->U{Q)+X 



Op 



dQ^dQldy K^tiQ'^ 



daAB^'y^(Q)+x daAB^-ii{Q)+x 



dQj^ dy 



dQj< dy 



(3) 



where the factor 2 is a consequence that the virtual photon has two transverse polarization 
states 0. 

When Qt and are both large, the differential cross section for producing a virtual 
photon of polarization A = T or L in Eqs. (0) and (0) can be factorized as 



da 



AB^^l{Q)+X 

dQ\ dy 



j dXiCpa/Ai^i, H) J dX2(l)b/Bix2,IJ^) 

a,b 



dan 



dQj. dy 



[xi,x2,Q,QT,y;fJ') (4) 



where J2a b over all parton flavors, the n represents both renormalization and factoriza- 
tion scales, and (pa/ a and (pb/B are normal parton distributions. The short-distance parts, 
daab^'yl{Q)x I dQ%dy in Eq. (|) are perturbatively calculable in terms of a power series in 
as ll^. When Qt ^ Q, the perturbatively calculated short- distance parts receive a large 
logarithm \\i{Q\/Q'^) for every power of beyond the leading order. Therefore, we need to 
resum such logarithms to all orders in in order to derive a reliable cross section for the 
virtual photon production at large Qt- 

It was shown [|lT| that the large logarithms {Q"^ / Q'^) come from the partonic sub- 
processes in which the virtual photon is produced from the decay of a parton that itself 
was produced at a distance scale ~ I/Qt- Because the parent parton can radiate massless 
partons before it radiates the virtual photon, there is one power of large \\i{Q\/Q'^) for an 
additional power of Ug- Since all powers of the large logarithms are from the final-state 
bremsstrahlung of a parton produced at a short- distance, such logarithms can be resumed 
into a fragmentation function for the parton to fragment into the virtual photon. Berger, 
Qiu, and Zhang [|l^] derived a modified factorization formula for calculating Drell-Yan cross 
section at Qt > Q, which includes the resummation of the large logarithms. 



da AB^'yl(Q)+x _ ^ r dz 

a,b,c 



dQT dy 



J2 J ^ J dXi(j)a/AiXi, H) J (ix2 06/B(x2,/i) 



X 



daab^cx 
dpl^ dy 



{xi,X2, Pc = Q/z;iJ, 



D^^^*{z,fi;Q) 



da^ 

+ I dXl(f)a/A{xi,fJ^) I dX2<f)b/B{x2,f^) J^^"^^ {Xl, X2, Q , Qt, V] f^) (5) 



a,b 



dQT dy 



where a, b, and c run over all parton flavors, Q'^ is the on-shell part of momentum Q^, and /i 
represents the renormalization, factorization, and fragmentation scales. The first term on the 
right-hand-side of Eq. represents the resummed contributions, which include all powers of 
the large logarithms; and the second term is free of the large logarithms and carries the same 
physical meaning as the F-term in the resummation formula for Qt <^ Q ||12|- In E<1- (H); the 



Dc--,-y* are the fragmentation functions for the partons of flavor c to fragment into a virtual 
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photon of invariant mass Q and polarization A, and their operator definitions are given in 
Ref. The resummation of the large \n{Q'^/Q'^) to all orders in can be achieved by 
solving the evolution (or renormalization group) equations for these fragmentation functions 
1^. In Eq. d^), the partonic part da ab^cx I dpl^dy are perturbatively calculable and represent 
the production of partons of fiavor c at a distance scale ~ 1/Pct ~ [0]- Since the 



F-term is defined to be the perturbative difference between Drell-Yan cross section and the 
resummed part of the cross section, the corresponding partonic parts at the order of a" are 



given by ||TT 



Y(n) 



= a 



in) 

ab^-rl {Q)X 



-E 



(m) 

ab^cX 



m=2 



in—m) 



(6) 



where ® represents the convolution over z as defined in Eq. and D*^"-' are perturbatively 
calculated virtual photon fragmentation functions at order of a". 

The modified factorization formula in Eq. (|^) should be valid for all Qt > Q- When 
Qt ^ Qi the resummed part dominates the cross sections, and the F-term represents a small 
correction. When Qt ~ Q-, the logarithms are small and the cross sections are dominated 
by the F-term. All short-distance parts in Eq. (|^) are evaluated at the same short- distance 
scale ~ 1/ Qt and free of the large logarithms. All leading logarithmic contributions from 
a distance scale of l/Qr to 1/Q are resummed into the fragmentation functions. With the 
modified factorization formula in Eq. , we are able to calculate the polarization parameter 
ody in Eq- (iD reliably in QCD perturbation theory by calculating the short-distance parts 
order-by-order in as- 

From the modified factorization formula in Eq. (^, we obtain the leading order contri- 
butions to the virtual photon cross sections by calculating the lowest order contributions to 
both dab-^c and a^^^.. From Eq. (^, the lowest order contributions to the short-distance 
d&ab^cx/dpl dy are independent of the virtual photon's polarization, and are given by the 



lowest order two-to-two partonic tree diagrams at 0{al 



da, 



ab—fcd 



dp^ 



CT 



dy 



Xl,X2, Pc = Q/z;fi) = TT 



Xi X2 S 



->cd 



6{s + i + u) 



(7) 



where g is the strong coupling constant and S = {Pa + PbY is the collision energy. In 
Eq. (^, s, t, and u are the parton level Mandelstam variables. The lowest order two-to-two 



matrix element squares in Eq. (|^) are given in Ref. []T3 

The lowest order contributions to the partonic cross sections aab^-^i in Eq. (H) are at order 
of 0{aemOis) from the subprocesses: qq l*{Q)g and qg l*{Q)q- Since the subtraction 
term in Eq. (|^) starts at order of 0{aemOi1), we have the leading order contributions to the 
F-term as 

.SY-LO) ^^(LO) 



da 



dQ"^ dy 



TT 



ab-^'yl(Q) 

dQ"^ dy 

Xi X2 S 



-M 



eg 



ab- 



6{s + i + u- Q'^ 



(8) 



where e is the electron charge and s, t, and u are again the parton level Mandelstam variables. 
The matrix element squares in Eq. (H) for producing a polarized virtual photon depend on 
the photon's polarization vector e^(Q). 
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The choice of the polarization vectors are not unique But, it should be consistent 
with the definition of ody in Eqs. (|l|) and (|]). The angle 9 in Eq. (|l|) depends on the frame 
choice and its 2;-axis |p. As a result, the extracted value of ody as well as the polarization 
vectors used to define the polarized virtual photon cross sections in Eq. are also sensitive 
to the frame choice and its z-axis. For the consistency of the QCD factorization formula in 
Eq. (^, it is very important to use the same polarization vectors to calculate both the virtual 
photon fragmentation functions and the F-term. Since we are interested in the region where 
Qt ^ Qi we will not choose the Collins-Soper frame |T^, which offers special advantage 
when Qt <C Q. Instead, as shown in Fig. |^, we choose the z-axis along the direction of 
the virtual photon's momentum Q and angle 9 to be the angle between the Q and Z, the 
direction of one of the decay leptons. Our choice of the 2;-axis corresponds to the "S-helicity" 
frame defined in Ref. With our choice of the helicity frame and angle 9, ody = 
corresponds to the virtual photon in a purely transverse (longitudinal) polarization state. 

To minimize the frame dependence, we present our results in terms of the covariant 
hadronic tensors W[ab]nv and polarization vector e^((5), which are defined as 



eg ^ 



iQ)exiQ) 



(9) 



For quark-antiquark annihilation subprocess: q{pi)q{p2) l*{Q)g-, we obtain the hadronic 
tensor, 



w 



i- 



1 

iu 



with quark's fractional charge e^; and for the "Compton" subprocess: g{pi)q{p2) 
we have 




(10) 



w 



[gq] 



s{-u) 



P2 ■ Q 



4Q^ 



Q 



u\ 



(11) 



With our choice of the helicity frame, we have the polarization vector for a longitudinally 
polarized virtual photon as [3 
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'XQ) = - ^ 



{P-QY-Q^S 



P-Q 



(12) 



where = P4 + Pg. The polarization vector e^Q) in Eq. (|I2|) reduces to the unit vector 
z'^ in the photon's rest frame |^. By contracting the hadronic tensors in Eqs. ( p!0D and 
(|rT|) with the polarization tensor ^*l^{Q)^'l{Q) (or —g^'^), one can easily obtain the matrix 
element squares for producing longitudinally polarized (or unpolarized) virtual photons. 
Because of the ambiguities in choosing the x- or y-axis, the functional forms of the transverse 
polarization vectors e^.. (Q) with i = 1,2 are not unique. One can derive the matrix element 
squares for producing a transversely polarized virtual photon by using the polarization tensor 
lJ2i=i,2^T!{Q)^TiiQ) 0! or by taking the one half of the difference in the matrix element 
squares for unpolarized and longitudinally polarized virtual photons. 

For deriving our numerical results, we use CTEQ5M parton distributions |]T3[ and the 



virtual photon fragmentation functions from Ref. |Q. We set the renormalization, factoriza- 



tion, and fragmentation scales equal to /i = /tyQ^ + Qj- with k a constant of order one. 

With the large Qt setting up the scale of hard collision, the QCD factorization formula 
in Eq. (^) should be valid for as small as a few GeV^. Since we are interested in the 
kinematics similar to J/ip production at high Q^, we choose Q = 2 and 5 GeV in following 
plots. 

In Fig. 1^, we plot the polarization parameter ody in the proton-antiproton collisions as 
a function of Qt at = 2 TeV (new Tevatron energy). The dashed lines represent the 
aoY with only the leading order F-term derived from Eqs. ([lOD and (pT]). The solid lines 
are equal to the ody with both the F-term and the resummed contributions. In Fig. ^, we 
plot the same polarization parameter any in the proton-proton collisions at \/S = 500 GeV 
(the RHIC energy). The constant k = 1 in both Figs. |^ and ^. By varying the k from 1/2 



to 2, we find that the cross sections are not very sensitive to the choice of k [|T1|, and the 
numerical values of the ody are relatively stable. 

From Figs. ^ and ^, we see that the virtual photons in Drell-Yan massive lepton-pair 
production are likely to be transversely polarized at high Qt- The resummed contributions 
provide significant corrections to the numerical values of the polarization parameter a^Y, 
and they are more important at the low Q^. As pointed out in Ref. @], the longitudinally 
polarized virtual photon fragmentation functions dominate the threshold region. Conse- 
quently, we expect that the resummed contributions to the cross sections of longitudinally 
polarized virtual photons are more important in low Qt than high Qt region, which are 
clearly evident in the difference between the solid and dashed lines. Since the perturbatively 
calculable partonic hard parts in the modified QCD factorization formula in Eq. (^ are 
evaluated at a single hard scale Qt without large logarithms, high order corrections should 
not alter the polarization parameter a^y dramatically fill]. We therefore conclude that the 
virtual photons of invariant mass Q in Drell-Yan massive lepton-pair production are likely 
to be transversely polarized when Qt 3> Q. 

When Qt ^ Mj/^, hadronic J/ip and the virtual photon production have a lot in com- 
mon. The key difference is the difference between their respective fragmentation functions, 
as shown in Fig. |1]. As argued in Ref. 0], the virtual photon fragmentation functions are 
completely perturbative, while the parton to J/ ijj fragmentation functions involve final-state 
radiations and soft interactions As shown in Fig. |l|, the fragmentation functions from a 
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parton d to a physical J/V' can be approximated by the fragmentation functions to a virtual 
gluon of invariant mass Q, which immediately decay into a cc-pair, convoluted with a transi- 
tion from the cc-pair to a physical J/V' 0. The first part of the fragmentation functions for 
a parton to a virtual gluon should be perturbatively calculable. Since a gluon can directly 
couple to a gluon, the evolution equations for the virtual gluon fragmentation functions can 
have a nonvanish leading order inhomogeneous evolution kernel for a gluon to a virtual gluon 
of mass Q. Consequently, the quark-to-virtual-gluon and gluon-to-virtual-gluon fragmenta- 
tion functions should be equally important, while the gluon-to- virtual-photon fragmentation 
functions are much smaller than quark-to- virtual-photon fragmentation functions How- 
ever, this difference should not have too much effect on the ratio of the transverse and 
longitudinal polarizations. We then expect that at high Qt, the virtual gluon, so as the 
cc-pair immediately produced from the decay of the virtual gluon, are more likely to be 
transversely polarized. 

The produced charm and anticharm quark pair of invariant mass Q can in principle 
radiate gluons and have soft interactions with other partons in the collisions. However, due to 
the heavy quark mass, such final-state interactions during the transition from the produced 
cc-pair to a physical ^/ip meson are not expected to significantly change the polarization 
0. If the formation from the cc-pair of invariant mass Q to a physical ^/if) meson does not 
change the polarization, one can expect the polarization of the ^/ip mesons produced at high 
Qt to be similar to the polarization of the virtual photon in Drell-Yan massive lepton-pair 
production at the same kinematics. For example, the non-relativistic QCD (NRQCD) model 
oi i /ip production precisely predicts the i /ifj mesons to be transversely polarized at large 
Qt 0, which is not consistent with recent Fermilab data 

Since i /ip production at high Qt involves both perturbative and nonperturbative physics, 
we have to find answers to the following two questions in order to resolve this inconsistency. 
The first question is about the reliability of the perturbative calculations for the relevant 
kinematics. The second question is about the separation between perturbative and nonper- 
turbative physics. These two questions can be summarized into one: if there is a reliable 
QCD factorization for hadronic ^/ip production @]. The measurements of the virtual pho- 
ton polarization in Drell-Yan massive lepton-pair production can help us to answer the first 
question. Since the QCD factorization for Drell-Yan massive lepton-pair production is ex- 
pected to be valid when Qt ^ Q [0], the measurements of the polarization parameter 
ody provide excellent tests of the reliability of perturbative calculations for the kinematic 
region similar io ^/ip production. Although the measurements of the a-^Y cannot test the 
final-state effect of J/V' production or answer the second question, the measured difference 
between ody and aj/^ can help us to isolate the role of the final-state effect, and to narrow 
the questions about the i /tp production. 

This work was supported in part by the U.S. Department of Energy under Grant No. 
DE-FG02-87ER40731. 
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FIGURES 





FIG. 1. Sketch for Drell-Yan massive lepton-pair and i/ip production via parton fragmentation. 
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FIG. 2. Definition of the angle 9 in Eq. (|) for Drell-Yan massive lepton-pair production. 
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FIG. 3. ody as a function of Qt in proton-antiproton collisions at ^/S = 2 TeV and Q = 2 GeV 
(a) and Q = 5 GeV (b). Dashed lines correspond to the Y-term only, and solid lines include both 
the y-term and the resummed contributions. 
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FIG. 4. q;dy as a function of in proton-proton collisions at ^fS = 500 GeV and Q = 2 GeV 
(a) and Q = 5 GeV (b). Dashed lines correspond to the y-term only, and solid lines include both 
the y-term and the resummed contributions. 
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